#for i in [2,4,8,16,32,64,128,256,512,1024,2048]
for i in [4096]
    global nels = i
    include("main1.jl")
    println(neq, "     ", res5[1, 2:end])
end

res5 = post1(20)
#=
plt.plot(res5[:,1], res5[:,2], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,3], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,4], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,5], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,6], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,7], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,3] .+ res5[:,6], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,4] + res5[:,7], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()

np.savetxt("results\\IntForceData.txt", res5)
=#
